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We introduce discrete systems in the form of straight (infinite) and ring-shaped chains, with 
two symmetrically placed nonlinear sites. The systems can be implemented in nonlinear optics 
(as waveguiding arrays) and BEC (by means of an optical lattice). A full set of exact analytical 
, solutions for symmetric, asymmetric, and antisymmetric localized modes is found, and their stability 

■ is investigated in a numerical form. The symmetry-breaking bifurcation (SBB), through which the 

asymmetric modes emerge from the symmetric ones, is found to be of the subcritical type. It is 
\ transformed into a supercritical bifurcation if the nonlinearity is localized in relatively broad domains 

around two central sites, and also in the ring of a small size, i.e., in efi'ectively nonlocal settings. 
The family of antisymmetric modes does not undergo bifurcations, and features both stable and 
I unstable portions. The evolution of unstable localized modes is investigated by means of direct 

simulations. In particular, unstable asymmetric states, which exist in the case of the subcritical 
bifurcation, give rise to breathers oscillating between the nonlinear sites, thus restoring an efi'ective 
^0 ' dynamical symmetry between them. 
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I. INTRODUCTION 

> ■ 

\ Spontaneous symmetry breaking is a fundamental effect caused by the interplay of nonlinearity With linear potentials 
■ featuring basic symmetries, such as double- well structures. In particular, while it is commonly known that the ground 
state in one-dimensional quantum mechanics follows the symmetry of the underlying double- well potential the 
self-attractive nonlinearity added to the respective Schrodinger equation [which transforms it into the Gross-Pitacvskii 
Jif^ ' equation for a Bose- Einstein condensate (BEC) of interacting atoms loaded into the double- well potential 0], or into 
the nonlinear Schrodinger equation in optical counterparts of the system 0] breaks the symmetry of the ground 
state, replacing it by a new asymmetric state minimizingthc system's energy, provided that the strength of the self- 
attraction exceeds a certain critical value (see, e.g., Ref. [1] for the general consideration, and Ref. 0] for the analysis 
of the symmetry-breaking self-trapping in BEC). The spontaneous symmetry breaking in double- well potentials was 
realized experimentally in BEC [6| and in nonlinear optics, where the symmetry breaking was observed in a setup 
' based on a photorefractive material Q. 

The symmetry-breaking bifurcation (SBB), which destabilizes the symmetric ground state and gives rise to an 
asymmetric one in the nonlinear system, was originally discovered in a discrete model of self-trapping [8| . In nonlinear 
optics, a similar SBB was predicted in Ref. 0] for continuous-wave (spatially uniform) states in the model of dual- 
core fibers. For solitons in dual-core systems, this bifurcation was studied in detail in Ref. Subsequently, the 
SBB was studied for gap solitons in dual-core fiber Bragg gratings with the same Kerr (cubic) nonlinearity as in the 
ordinary fibers [Tl] | . In another physical setting which also features the cubic nonlinearity, the SBB was predicted for 
matter- wave solitons in the self- attractive BEC loaded into a dual-trough potential trap jl^-[l5|. 

The self-focusing cubic nonlinearity gives rise to the soliton bifurcations of the subcritical (alias backward) type, 
in which the branches of asymmetric modes emerge as unstable ones, going backward and getting stabilized after 
switching their direction forward at turning points [l6j . On the other hand, the combination of the self-focusing 
nonlinearity with a periodic potential acting in the free direction (perpendicular to the direction of the action of the 
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double-well potential) changes the character of the bifurcation from subcritical to supercritical, with the asymmetric 
branches emerging as stable ones and immediately going in the forward direction [12, US] • The SBB in the model of 
the dual-core fiber Bragg grating is of the forward type too [ll|. The sub- and supercritical (alias backward/forward) 
SBBs may be regarded as examples of phase transitions of the first and second kinds, respectively [l6} . 

A physically interesting alternative to the linear double- well potential is the setting with an effective pseudopotential 
[l3| induced by the double^eak spatial modulation of the local nonlinearity coefficient, which may be implemented in 
optics and BEC alike [3 [13 • The ultimate form of such a setting is the one with the nonlinearity concentrated at two 
points, in the form of a symmetric pair of delta- functions or narrow Gaussians [20l. [2l|. as well as a two-dimensional 
counterpart of the system, in the form of two parallel troughs in which the nonlinearity is applied [l^ (the underlying 
model, with the self-attractive nonlinearity represented by a single delta-function, was introduced much earlier in Ref. 
[1^). The SBB of solitons in the symmetric double- well nonlinear pseudopotentials was studied recently, featuring 
the subcritical type of the symmetry breaking [20l - [2^ . 

The spontaneous symmetry breaking was also analyzed for solitons in dual-core discrete systems, with the uniform 
coupling between two parallel chains [l^l, or with the coupling established at a single site [IBl- In the former case, 
the SBB is subcritical, while in the latter case it is supercritical. 

The objective of the present work is to consider the SBB in one-dimensional discrete lattices with the nonlinearity 
tightly concentrated at two symmetric sites, or in narrow regions around them. This is a straightforward discrete 
counterpart of the double- well nonlinear pseudopotential [lo, [HI, which offers a simple testbed for the study of 
symmetry-breaking effects in discrete media. Physically, the linear chain with two nonlinear sites can be readily 
implemented in optics, by embedding two nonlinear cores into an arrayed linear waveguide, and in BEC, by means 
of the Feshbach-resonance technique applied locally to the condensate trapped in a deep optical lattice 0- An 
essential advantage of the system is that, as we demonstrate below, it admits a fully analytical solution (for the 
infinite chain), with an arbitrary separation between the two symmetric nonlinear sites, while the stabilit y o f the 
exact solutions may be efficiently predicted following general principles of the elementary bifurcation theory [l6| and 
the Vakhitov-Kolokolov (VK) criterion [26j . 

The paper is organized as follows. The model is introduced in Section 2, which is followed by producing exact 
solutions for symmetric, asymmetric, and antisymmetric localized modes in Section 3. Numerical results, obtained 
for finite lattices, are reported in Section 4 (in particular, it is demonstrated that the SBB is subcritical in the system 
with two nonlinear sites). In Section 5, we consider an essentially different version of the system, in the form of a 
ring-shaped chain, with the two nonlinear sites placed at diametrically opposite points, in which case an analytical 
solution for the SBB is available too. The paper is concluded by Section 6. 



II. THE MODEL 



According to what is said above, the model is based on the linear discrete Schrodinger equation with two nonlinear 
sites embedded into it: 

.du^ + (1/2) (u„+i + U„_i - 2Un) + {5n,0 + Sn.l) \u„f U„ = 0, (1) 

at 

where / is the integer distance between the two nonlinear sites, and 6n^m is the Kronecker's symbol. The evolutional 
variable, t, is time in the application to BEC, or propagation distance in the case of an array of optical waveguides. In 
the former case, the nonlinearity at two sites can be induced by focusing laser beams, that may enhance the nonlinearity 
through the Feshbach resonance [l^l, at two particular droplets of the condensate trapped in a deep optical lattice. 
In the photonic realization of the setting, strong nonlinearity in two particular cores in the waveguiding array can be 
readily imposed by doping them with resonant atoms (see, e.g., Ref. [l^l), which does not affect the linear coupling 
of these sites to adjacent ones, as implied in Eq. ([T}. The same mechanisms can be used for inducing the local 
nonlinearity in circular configurations considered below in Section 5. 

Stationary solutions to Eq. ([T]) are sought for as u„(t) — e~*"*C/„, where real stationary field Un obeys equation 

LuUn + (1/2) (Un+l + Un-1 - 2C/„) + {5,,fl + 5,-,,i) uf, = 0. (2) 

Along with this model, we will also consider its version with a smoothed form of the discrete delta- function, viz., 

Sn,no exp ^- (n - nof / A^^ , (3) 
where no = or ^, and A is the smoothing width. 
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LoU + -— + 26{x)U'^ ^Q, (7) 



A preliminary remark is that Eq. ([2]) with I = corresponds to the single nonlinear site with the double strength, 
placed at n = 0. In that case, an obvious exact solution is 

(C^n)singic = ^e-'^l"! , ^ (1/2) sinh «, (4) 

where k is connected to uj by the dispersion relation for evanescent waves in the linear lattice: 

w= -2sinh2(K/2). (5) 

Taking relation ([5]) into regard, the norm of solution (j4]) is 

N = C/^ ^ (1/2) cosh K= (l-u;)/2. (6) 

n—~QO 

Note that this expression for N satisfies the VK stability criterion, dN/duj < 0, hence solutions Q may be stable 
[2^. On the other hand, the continuous counterpart of Eq. ([2]) with the single nonlinear site is 

2 

where S{x) is the delta- function (this continuous equation was first introduced in Rcf. [23|). It has an obvious localized 
solution, 

U{x) = exp (-^/^|a;|) , (8) 

whose norm is degenerate (it does not depend on w): N = J^°° U^{x)dx = 1/2. Being formally neutrally stable in 
terms of the VK criterion, all solutions are actually unstable f2ll| . The degeneracy and instability of solutions ^ 
resemble the classical properties of the Townes' solitons (30| . 

III. EXACT SOLUTIONS FOR THE INFINITE LATTICE 

A. General analysis 

Symmetric, antisymmetric, and asymmetric solutions to Eq. ^ can be found in an exact form, following the 
pattern of the exact solutions for the continuous counterpart of Eq. ^ , which was considered in Ref. [l^l : 

L^U+l^ + [Six)+5{x~l)]U^ ^0, (9) 

cf. Eq. Exact symmetric and asymmetric solutions to the discrete equation ^ are sought in the following form: 

( Ae™, at n < 0, 

Un = I B cosh (k (n - no)) , at < n < (10) 
[ Ce-''("-'), a.tn>l. 

Coordinate hq determines the location of the center of the intermediate part of the solution. Note that uq does not 
need to be an integer number. Symmetric and asymmetric modes correspond, respectively, to hq = 1/2 and uq ^ 1/2. 

Ansatz automatically satisfies the linear discrete Schrodinger equation. There remains to check Eq. ([2]) at the 
nonlinear sites (of course, imposing the condition of the continuity of the solution at these sites). With regard to Eq. 
([5|), the continuity condition yields relations between the amplitudes: 

A C 
cosh(Krio) cosh {kI — kuq) ' ^ 

and the equation at the nonlinear sites amounts to the following relations: 

A' = , = ""!\" (12) 

After some algebra, the condition that Eqs. (fTTj) and ([T2|) yield the same expression for B produces an equation for 
no: 

x^ + {3- L) a;^ + (1 - SL^^) x - L'^ = 0, (13) 
a; = e-2™«, L = e^^K (14) 
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B. Symmetric modes 



Two roots of quartic equation p3| are x = ±L-i/2. The negative one is unphysical, while the positive root 
corresponds, according to Eq. (jl4p . to uq = 1/2, which represents the symmetric mode. The amphtudes of the 
symmetric solution, as given by Eqs. (ITT|) and are 



A' 



sinh K 



4e sinh k 



(1 



(15a) 
(15b) 



It is relevant to note that, in the case of Z = 0, Eqs. (|15p yield ~ ~ ^ (1/2) sinhK, which naturally coincides 
with as given by Eq. 



C. Asymmetric modes 

Two other roots of Eq. p3p represent a pair of asymmetric modes: 



1 



x± 



(L-3)±(L-1) 



L 



(16) 



Taking into regard Eq. (|14p . it is easy to check that roots ([16]) correspond to two values ('t.o)± which are located 
mutually symmetrically around the center: a;+x_ = i^^, i.e., (f^o)-|- + ('^o)_ = ^- Further, roots (|16p are physical if 
they are real and positive, which means L > 4. Thus, the symmetry-breaking bifurcation (SBB), i.e., the appearance 
of the asymmetric modes from the symmetric one [see Fig. [IJa)], happens, with the increase of distance I between the 
nonlinear sites (at fixed k, i.e., fixed w), at the critical point, KcJ — In 2, or, in terms of the frequency of the localized 
mode, at 



-2 sinh^ ,' 

21 



(17) 



In other words, for fixed I, the SBB occurs with the increase of |aj| at point (|17p . the asymmetric modes existing at 

1^1 > kcl-l- 

The asymmetric mode keeps the double-peak shape if the minimum point in expression ()10p . hq, is a really existing 
minimum, rather than a virtual one, i.e., it stays in interval < no < /. Further, this means that x± must fall into 
the following interval: L^^ < a;± < 1. It is easy to check, making use of Eq. (|16p . that the latter condition places L 
and / into narrow intervals of their values, namely, 4 < L < 2 + \/5 ~ 4.24, and, accordingly. 



In 2 w 0.69 < < (1/2) In (2 + Vd) ^ 0.72, 



(18) 



cf. Eq. (dH). At kI > (1/2) In (2 + \/5) , i.e., cj < cj+ = -2 sinh^ (Al) ^ In (2 + [in other words, at |a;| > 

cf. Eq. pT]) ]. no leaves the region of < no < /, hence the asymmetric mode becomes single-peaked ("strongly 
asymmetric"). The transition from the double-peak asymmetric mode to the single-peak shape is illustrated by Fig. 

Elb). 

The final expressions for the amplitudes of the asymmetric mode arc obtained by the substitution of expressions 
(ITB and ([H]) into Eqs. ^ and (|T^ : 




26**' sinh K 



- 3) ± (e^-^' - 1) Ve2«' - 4 



sinh K 



(19a) 
(19b) 
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where signs ± pertain to the two conjugate asymmetric modes. Note that, at the SBB point, e = 2, expressions 
(|19p coincide with their counterparts (jlSp obtained above for the symmetric mode: at this point, the amplitudes are 
= C^ = (2/3)sinhK and = (16/27) sinh k. 
As a natural measure of their asymmetry, we will use 

© ^ A2-fc^ = ±Vl-4e-2«'- (20) 

Obviously, = at the SBB point, e"**' = 1/2, and expression ([20|) makes sense past the bifurcation point, i.e., 
at e^"' < 1/2. In particular, at the above-mentioned point of the transition from the double- to single-peak profile, 
kI = (1/2) In (2 + \/5), the asymmetry is still relatively smah, |6| = - 2 w 0.24. 



D. Antisymmetric modes 

Antisymmetric solutions, and their counterparts with broken antisymmetry (if any) are looked for by dint of the 
following ansatz, cf. Eq. (fTO)) : 



Ae'=", at n < 0, 

Un ~ i B sinh {k {n — no)) , at < 7i < /, 



(21) 



Looking for a bifurcation occurring to the antisymmetric mode, one arrives at the respective roots for x [with x defined 
as in Eq. (fT4)) ]: 



1 



(L-3)±(L-1) 



L 



(22) 



which differs from Eq. (jl6p by the opposite sign in front of the first term in the square brackets. This difference makes 
both roots (|22p negative (unphysical), hence the antisymmetric solutions do not undergo the bifurcation, similar to 
the situation known in the continuous model (poj). 

The mode with the unbroken antisymmetry is given by the exact solution in the form of ansatz ([21]) with amplitudes 



A = -C = 



sinh Hi 



B 



1 - e^'^^ ' 
2\/sinhKe~'^'/^ 
(l-e-^^)^^^ ■ 



(23a) 
(23b) 



Finally, it is natural to call the antisymmetric modes, corresponding to even and odd I, on-site and inter-site ones, 
respectively, as they have the zero point, hq = 1/2, either coinciding or not with the (central) site of the lattice. 



IV. NUMERICAL RESULTS FOR A FINITE LATTICE 



The objective of the numerical solution is to obtain solutions of stationary equation ^ for a finite lattice, compare 
them to the exact solutions found above for the infinite lattice and identify stability of the symmetric, asymmetric 
and antisymmetric solutions. In addition, the numerical calculations help to identify the type of the SBB (sub- or 
supercritical), as the analytical solution produces a very cumbersome result, in this respect. Numerical stationary 
solutions were constructed in the lattice of 71 sites. Highly asymmetric modes were obtained by means of the 
continuation in w, using the Newton's iteration procedure that started from the analytical asymmetric solution given 
by Eqs. pU]) and p^ . For all numerical solutions, we analyzed the linear stability by solving the corresponding linear 
eigenvalue problem. The results were checked by direct simulations of the underlying equation ([1]). 



A. Symmetric modes 



In Fig. [21 families of numerically generated symmetric solutions are shown for different distances (/) between the two 
nonlinear sites. As said above, the stability of the numerically found solutions was identified through the calculation 
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FIG. 1: (Color online) (a) The bifurcation diagram in the plane of the intrinsic frequency of the localized modes {to) and their 
total norm (A'^), as produced by the analytical solution for the infinite lattice with two symmetric nonlinear sites, separated 
by distance I = 4. The black thick, red thin and blue dotted lines correspond to the symmetric (|15|l . antisymmetric (|23|) 
and asymmetric (|19|l solutions, respectively. In the insert, the zoom of the bifurcation of the branch of asymmetric solutions 
(the blue dotted line) is shown. The gray area determines the region of the existence of the two-peaked asymmetric solution, 
limited by critical frequencies < uj < ujc^. In panels (b)-(d), examples of analytically found symmetric, antisymmetric and 
asymmetric modes at a; = —0.016 and I = 4 are shown. 

of the corresponding eigenvalues of small perturbations. The borders between stable and unstable segments of the 
solution branches correspond to the SBB which destabilizes the symmetric solutions. The numerical solutions are 
practically indistinguishable from their analytical counterparts (|10p . (|15p . therefore the analytically found curves are 
not plotted separately (strictly speaking, the numerical results cannot be identical to the analytical ones, as the 
numerical computations were performed for the finite chain, while the analytical findings pertain to the infinite one) . 

To check the predictions of the linear stability analysis, we simulated the underlying equation ([1]), with initial 
profiles for symmetric modes taken in regions where these solutions are expected to have different stability. Typical 
examples, presented in right panel of Fig. [21 show that the unstable symmetric stationary solution transforms into a 
pulsating single-peak mode, which breaks the symmetry, getting spontaneously localized on one of the nonlinear sites. 

B. Asymmetric modes 

The numerical results for the asymmetric modes are displayed in Fig. [3ja). In particular, unstable parts of the 
families of asymmetric solutions are those which are related to the SBB of the subcritical type [see Fig. IDJb) below]. 
The SBB displayed by the numerical results closely follows the analytical solution for the infinite lattice. In particular, 
generic numerically found profiles of the double-peaked (above w^) and single-peak (below w^) asymmetric modes 
are compared with their analytical counterparts in Fig. [3l^b), the difference between the numerical and analytical 
solutions being < 1%. 

In Fig. m the numerical results obtained for the asymmetric solutions are collected in the form of plots showing 
the asymmetry measure defined, as per Eq. (j20p . through the difference between squared amplitudes of the solution 
at the two nonlinear sites, as a function of frequency co and total norm N. The corresponding discrepancy between 
the numerical and analytical results is < 1%. In panel (a), the numerically found bifurcation points are compared to 
those predicted analytically by Eq. (jl7p (they correspond to values of a;~ indicated in the panel), which demonstrates 
a precise agreement. An important conclusion clearly suggested by Fig. Sl^b) is that the bifurcation has the subcritical 
character, with the asymmetric branches being unstable exactly between the SBB and turning points, as might be 
expected. 

To check the predictions of the linear-stability analysis, we simulated the underlying equation ([!]) with initial profiles 
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FIG. 2; (Color online) Left panel: Families of numerically generated symmetric solutions (their analytical counterparts are 
indistinguishable from the numerical ones). The corresponding values of the distance between the nonlinear sites, I, are 
indicated. The solid and dashed lines depict, respectively, stable and unstable portions of the solution families. Right panel: 
The evolution of stable (a) and unstable (b) symmetric modes with / = 4, taken at uj = —0.014 and uj — —0.02, respectively. 
Here and in similar plots below, the spatiotemporal evolution is displayed by means of the density contour plots. Small 
perturbations (~ 1%) were added to initiate the evolution. Dashed lines show the position of nonlinear sites. 




FIG. 3: (Color online) (a): The bold and thin lines show, respectively, the families of asymmetric and symmetric solutions. As 
in Fig. [2j solid and dashed lines corresponds to stable and unstable portions of the solution families. In fact, the lines for the 
symmetric solutions are identical to those displayed in Fig. [21 (b): An example of the transition from the double-peak shape 
for ij > cj^ to the single-peak one for lj < tj^, obtained from numerical (circles and squares) and analytical (red dashed and 
black solid lines) solutions. 

for asymmetric modes taken in the regions where these solutions are expected to be stable and unstable, respectively. 
Typical examples, presented in Fig. [SJ show that the unstable asymmetric stationary solution transforms itself, after a 
transient period, into a robust breather oscillating between two asymmetric configurations (in that sense, the breathers 
restores an effective dynamical symmetry). 



C. Antisymmetric modes 

Families of numerically found on-site (a) and inter-site (b) antisymmetric modes are shown in Fig. [Bl As in the 
previous cases the curves, they are practically indistinguishable from the analytically found counterparts. In contrast 
to the symmetric and asymmetric solutions, norm N of the antisymmetric ones is bounded by a minimum value (the 
existence threshold). It is also seen that each curve features stable and unstable portions, the border between which 
approaches the bottom of the curve with the increase the distance between the nonlinear sites, I. Panels (c) and (d) in 
Fig. ini display typical profiles of the on-site and inter-site antisymmetric modes, again showing very good agreement 
with the respective analytically predicted profiles. Note that, according to panels (a) and (b), the antisymmetric 
solution shown for ^ = 4 is stable, while the ones for Z = 1, 2, 3 are unstable. 

To check the linear-stability predictions for the antisymmetric modes, we took initial profiles of the antisymmetric 
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FIG. 4: (Color online) The family of the numerically found asymmetric solutions in the planes of {0,ijj) (a) and (0,A'') (b) 
for different distances between the nonlinear sites (l). As above, the stable and unstable portions of the solution families are 
depicted by solid and dashed lines, respectively. 
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FIG. 5: ((Color online) The evolution of unstable (a) and stable (b) asymmetric modes with I — 4, taken at u = —0.02 and 
uj — —0.05, respectively. Dashed lines show the position of nonlinear sites. 



solutions from Fig. [6l^c),(d) at w = —0.5, when the on-site/intcr-site modes with I = 1,2,3 are unstable, while the 
one with ^ = 4 is stable, and simulated Eq. ([1]) with small initial perturbations (~ 1%). The results are presented 
in Fig. [71 The unstable solutions spontaneously transform into single-peak modes localized on one of the nonlinear 
sites, which demonstrates decaying oscillations of the amplitude. On the other hand, the linearly stable antisymmetric 
mode, with ? = 4, is indeed robust in the direct simulations. 



D. Effects of the finite extension of the nonlinear region 



In the model with two finite nonlinearity domains, defined as per Eq. ([3]), results were obtained in the numerical 
form. They are displayed for the families of symmetric and asymmetric modes in Fig. |8l and for the antisymmetric 
ones in Fig. [9l In particular, Fig. [8] demonstrates that the increase of the width of the nonlinearity domain transforms 
the subcritical bifurcation into a supercritical one, thus completely stabilizing the asymmetric modes and making the 
value of N at the SBB point lower. It is also worthy to note that, as seen in Fig. ^ the broader nonlinearity gives 
rise to an additional stability window close to the existence threshold (minimum value of N). 
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FIG. 6: (Color online) Panels (a) and (b) display N{uj) curves for families of the antisymmetric solutions of on-site and inter-site 
types, respectively (i.e., even or odd distance I between the nonlinear sites). The solid and dashed portions of the curves refer, 
as usual, to stable and unstable solutions. Panels (c) and (d) display typical profiles of the antisymmetric on-site and inter-site 
modes, as obtained in the numerical and analytical forms aX ui = —0.5. 
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FIG. 7: (Color online) The evolution of unstable (a), (b), (c) and stable (d) antisymmetric modes. Initial profiles are taken 
from Fig. [6l;c),(d) with w = -0.5, for I = 1,2,3,4. 



V. THE RING MODEL 



If the infinite linear chain is replaced by a circle (ring) with the nonlinear sites placed at diametrically opposite 
points, Eq. ^ is replaced by the following one: 

U^Un + \ {Un+1 + f/„-l - 2C/„) + ((5„,„_ + Sn,n^) = 0, (24) 

where n = n± may be defined as the top and bottom points of a vertical diameter cutting the circle. In fact, it is 
more convenient to replace Eq. (j24p by a system of two identical equations for two semi-circles, left and right ones. 
In each equation, n takes values n_ < n < 71+ (n_ < and n+ > are assumed). In BEC, the circular chain can be 
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FIG. 8: (Color online) (a) Curves N{ij) of the symmetric and asymmetric solution families for I = 4 with A = 0.1 in Eq. Q 
(black thin lines, which are actually tantamount to their counterparts with the 5-like nonlinearity) , and with A = 1 (red thick 
lines), which correspond to a relatively broad Gaussian. In (b) and (c), the corresponding asymmetry ratio O is shown versus 
u and A'^. 
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FIG. 9: (Color online) The N{io) curves of antisymmetric solution families for I = 4 with A = 0.1 (black curves) and A = 1 
(red curves), cf. Fig. [S] 

realized as a combination of a toroidal trap and periodic potential created in it [3l] |. and in optics it corresponds to 
an array of waveguides created in a hollow cylindrical shell, or written in the form of the ring in a bulk sample [32| . 

In principle, one may expect two types of the symmetry breaking in this setting: between the top and bottom, 
which is a counterpart of what was considered above in the framework of the rectilinear chain, or between the left 
and right semi-circles. 

The linear solutions for the left (l) and right (r) semi-circles can be looked for as 

(t^n);,r = ^l,r COSh {k [u - Ui^r)) , (25) 

with K related to uj by Eq. ([S]), and some constants Ai_r and ni^r {ni and are not necessarily integer numbers). 
Then, the continuity conditions should be imposed at points n — 7i±, where the two semi-circles are linked into the 
entire circle: 

Ai cosh {n± — ni) = Ar cosh {n± — rir) ■ (26) 

After simple manipulations, one may eliminate the amplitudes from two equations (|26p . which leads to the following 
consistency condition: cosh (n (n_|_ — n_ — ri; + rir)) = cosh {n (n_|_ — n_ + n; — rir)) . It is obvious that the consistency 
condition can be met in the case of n+ ~ n_ , which is trivial (zero length of the ring) , or 

ni = rir = TT-Q. (27) 

Further, it then follows from Eqs. that a consequence of Eq. (P7)) is Ai = Ar = A, hence the symmetry breaking 
between the left and right circles is impossible. 

However, the top-bottom symmetry breaking is possible. Substituting ansatz (1251) with the left-right symmetry into 
Eq. (|24p at points n ~ n±, the result of a straightforward analysis is a system of two equations, corresponding to + 
and — : 



tanh {k {n± — no)) [l — tanh^ (k {n± — no))] — / sinh n, 



(28) 
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where uq is defined as per Eq. ((27|) . In equations ([28]) . and uq are considered as two unknowns, for given k and 
n±. As follows from Eqs. (pS)) . the symmetric solution, with tiq = (1/2) (n+ + n_), has amplitude 

-4symm (sinhK) tanh ((k/2) (n+ - [l - tanh^ (('«/2) (7i± - no))] . (29) 

An obvious corollary of Eqs. ([28|) is relation 

tanh (k (n+ — no)) [l — tanh^ (k (n^ — no))] — tanh (k (no — n_)) [l — tanh^ (k (no — n_))] = 0. (30) 

The SBB point can be found by setting no = (1/2) (n_|_ — n^) + Suq, with infinitesimal Shq^ and demanding that the 
coefBcient in front of Suq in the respective expansion of the left-hand side of Eq. ([50]) vanishes. The result is that this 
happens at point cosh((K/2) (n+ — n_)) = ■\/3/2, which is equivalent to 

[K(n+-n_)]„, = ln(2 + \/3) . (31) 

At this point, the amplitude is 

^2^. = (2/33/2^ sinhK. (32) 

The asymmetric state exists, for given k, i.e., given ui, for the length of the semi-circle, n+ — n_, which is larger than 
the one corresponding to Eq. pi[) . Alternatively, at given — n_, the asymmetric state exists for k exceeding the 
value defined by Eq. ([5T|) . 

In Fig. [To] we summarize the results related to the ring configuration. In panel (a), families of the symmetric and 
asymmetric solutions are displayed for three different lengths of the ring. We have checked the linear stability of the 
corresponding solutions through the numerical computation of the eigenvalues for small perturbations. The respective 
stable and unstable regions of the existence curves are shown by the solid and dashed lines, respectively. As in the 
case of the linear chain, at critical frequency Wcr calculated from Eq. (j3ip . the branches of the asymmetric solutions 
bifurcate from the branch of symmetric solutions. The amplitude of the symmetric solutions at the bifurcation point 
is displayed in Fig. [TUT b) for different lengths of the ring. 

Also for different lengths of the ring, we have analyzed the type of the corresponding SBB, calculating the asymmetry 
parameter Q as per Eq. (j20p . The dependences of & on frequency uj and norm are shown in Figs. llOf c) and (d). 
As one can see, the length of the ring plays a crucial role in the determination of the bifurcation type. In the small 
ring (e.g., of length 6), the bifurcation is supercritical, while the increase the length to > 8 changes it into a subcritical 
one. 

In Figs. [11] and I12[ the evolution of stable and unstable asymmetric and symmetric solutions is shown. The 
direct simulations corroborate the linear-stability analysis. As one can see, the solutions taken on stable parts of the 
existence curves are stable indeed. The unstable asymmetric solution [see Fig. HHb)] starts to oscillate between two 
asymmetric configurations, while the unstable symmetric solution in Fig. [T2jb) rapidly transforms into an oscillating 
single-peak mode. It is relevant to stress that, similar to the situation in the straight chain (cf. Fig. [5]), the instability 
of symmetric modes transforms them into asymmetric ones, while the instability of asymmetric modes leads to the 
emergence of the effective dynamical symmetry. 

Finally, it is relevant to note that, in the case of the circle of a finite length, it is also possible to consider values 
id > 0, i.e., imaginary wavenumbers, k = z|k|, see Eq. ([S]). This means that the exponentially decaying discrete waves 
are replaced by oscillatory ones. After a simple analysis, the respective version of Eq. (|28p can be derived in the 
following form: 

tan (|k| (n± - no)) [l + tan^ (|k| (n± - no))] = T^V [2 sin (|k|)] . (33) 

A straightforward analysis of Eq. (|33p demonstrates that, on the contrary to Eq. (|28p . it does not give rise to the 
symmetry breaking (the crucial difference is the opposite sign in front of tan^ in the square brackets). 

VI. CONCLUSION 

The objective of this work is to elaborate the simplest setting for the study of the spontaneous symmetry breaking 
in dynamical chains. For this purpose, we have introduced two one-dimensional discrete systems, in the form of the 
straight and ring-shaped linear chains, with two symmetrically inserted nonlinear sites. This pair of sites introduces the 
symmetry that may be spontaneously broken. The chains, both rectilinear and circular ones, can be realized in BEC 
(with the help of optical lattices) and in optics, in the form of waveguiding arrays. A full set of analytical solutions for 
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FIG. 10: (Color online) (a): The family of symmetric (thin) and asymmetric (thick) solutions for three different lengths of the 
ring (6, 8, and 10). The solid and dashed lines correspond to the stable and unstable solutions, respectively, (b): The amplitude 
of the symmetric solution at the bifurcation point, for different lengths of the ring (indicated by numbers at the points), (c), 
(d): The dependence of asymmetry O on the frequency (c) and norm A'^ (d), showing the change of the bifurcation type in the 
ring with the increase of the ring's size. 



symmetric, asymmetric, and antisymmetric localized states has been obtained in the explicit form, for both geometries. 
The stability of the stationary modes was investigated through the numerical computation of the eigenvalues for small 
perturbations. The SBB (symmetry-breaking bifurcation), which is responsible for the asymmetric modes emerging 
from the symmetric ones, is of the subcritical type in the straight lattice, and in the circular one of a sufficiently large 
size. The bifurcation becomes supercritical if the system is made effectively nonlocal, i.e., in the ring of a smaller 
size, or if the nonlincarity is spread over relatively broad areas around the two central sites in the straight chain. 
The antisymmetric modes are not subject to bifurcations, although they too may be both stable and unstable. The 
development of the instability (when it occurs) was tested with the help of direct simulations. It was found that 
the unstable stationary asymmetric states, which are a part of the subcritical bifurcation, spontaneously transform 
into breathers oscillating between the two nonlinear sites, that may be considered as the restoration of an effective 
dynamical symmetry. 
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FIG. 11: (Color online) The evolution of stable (a) and unstable (b) asymmetric solutions with initial profiles taken at uj 
and u} = —0.05, respectively. The length of the ring is 10 sites. Initial profiles were taken with small perturbations. 
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FIG. 12: (Color online) The evolution of stable (a) and unstable (b) symmetric solutions with initial profiles taken at u) - 
and uj = —0.03, respectively. The length of the ring is 10 sites. Initial profiles were taken with small perturbations. 
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